Light WIMPs in the Sun: Constraints from HeHoseismology 
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We calculate solar models including dark matter (DM) weakly-interacting massive particles 
(WIMPs) of mass 5-50 GeV and test these models against helioseismic constraints on sound speed, 
convection zone depth, convection zone helium abundance, and small separations of low-degree p- 
modes. Our main conclusion is that both direct detection experiments and particle accelerators 
may be complemented by using the Sun as a probe for WIMP DM particles in the 5-50 GeV mass 
range. The DM most sensitive to this probe has suppressed annihilations and a large spin-dependent 
elastic scattering cross section. For the WIMP cross-section parameters explored here, the lightest 
WIMP masses < 10 GeV are ruled out by constraints on core sound speed and low-degree frequency 
spacings. For WIMP masses 30-50 GeV, the changes to the solar structure are confined to the inner 
4% of the solar radius and so do not significantly affect the solar p-modes. Future helioseismology 
observations, most notably involving g-modes, and future solar neutrino experiments may be able 
to constrain the allowable DM parameter space in a mass range that is of current interest for direct 
detection. 



INTRODUCTION 

We explore the role of DM WIMPs in modifying the ther- 
mal gradient of the Sun. A similar study involving stan- 
dard WIMPS of mass 50 GeV or larger has been per- 
formed by Bottino et al. [I]. However, here we wish to 
consider the effects of low mass WIMPs, with masses as 
low as 5 GeV, with large trapped abundances within the 
Sun. 

To affect the Sun's thermal gradient, we need large 
elastic scattering rates. The solar sound speed can be 
affected, and helioseismology has been proposed as pro- 
viding a possible constraint on supersymmetric WIMPs 
[2^4]. For the range in masses in which we are inter- 
ested, the limits on the size of spin-independent WIMP 
elastic cross sections from CRESST [5 , XENON-10 [5], 
and XENON- 100 [7] are already quite stringent making 
it unlikely for helioseismology to provide further restric- 
tions. Moreover, recent COUPP results [5] (especially 
at masses > 10 GeV) and the results from PICASSO [9] 
similarly restrict spin-dependent interactions. However, 
these limits become weaker as the DM mass is decreased, 
especially for masses of around 5 GeVor less. 

Detailed solar models [1] show that a WIMP signal 
is only possible for cross sections in a limited range. 
For a DM particle mass of 50 GeV the relevant effective 
cross section for these signals was found to be of order 
10~^^ cm^. The limits on spin-dependent scattering from 
direct detection experiments, e.g., GOUPP [SUTU] restrict 
the spin-dependent elastic scattering cross section for a 
50 GeV DM particle to less than ~ 10"^'^ cm^. Moving to 
lighter masses alleviates these limits somewhat, and for 



around 5 GeV we are in the interesting region of around 
10~^^ cm^ [3]. In addition, given the astrophysical un- 
certainties that can affect these limits (e.g. [TT]) it is in- 
triguing to ask whether helioseismology can complement 
direct detection limits at lower masses. 

Recently, there has been interest in exploring the low 
DM mass regime as a possible way to consistently com- 
bine the results from the DAMA direct detection exper- 
iment [12] with those from others such as GDMS [13] 
and CoGeNT [H] (see e.g., [B]). Although we do not 
attempt to do the same here, we simply note that this 
mass regime is of great interest with upper limits on the 
spin-dependent elastic scattering cross section for low- 
mass DM, reaching ^ 10"'^^ cm^ ^16j for certain models 
and assumptions. 

Solar effects are most pronounced for DM with a sup- 
pressed annihilation cross section such that after capture 
by the Sun, the DM candidates do not annihilate quickly. 
A prominent example of this is asymmetric DM where 
annihilation is completely suppressed. 

In this paper, we outline a class of models for WIMPs 
that are capable of modifying the temperature profile in 
the core of the Sun, and illustrate their effects on he- 
lioseismology and neutrino fluxes. The accumulation of 
these WIMPs in the solar core results in significant en- 
ergy transfer to solar protons. We note that the effect is 
not large enough to account for discrepancies between ob- 
servations and helioseismology for models that also pre- 
dict the observed neutrino flux. However, given the cur- 
rent debate about the appropriate element abundances 
to be adopted in solar models [17], this effect may still 
play a role and should be included in the models. Indeed 



the recently revised solar abundances result in solar mod- 
els that cannot reproduce currently observed helioseismic 
data [181IT9]. 

The effects on the Sun of low- mass asymmetric WIMPs 
possessing large self-interactions have been considered in 
PP] , While in HO] the authors focus on WIMPs with 
spin-independent interactions, in this study we specifi- 
cally focus on WIMPs with purely axial interactions and 
consequently only spin-dependent elastic scattering. 

In the following sections, we outline the main features 
of DM models that can be potentially probed using solar 
properties and explore their effect on solar models. 



WIMP MODEL CHARACTERISTICS 

The most interesting DM mass region is the low mass 
region around 10 GeV and below. As stated in the intro- 
duction, direct detection constraints restrict the size of 
the DM elastic scattering cross section. In particular, the 
constraints on spin-independent elastic scattering for low 
masses by CRESST [S] are too stringent for any possible 
improvements from helioseismology. We will therefore fo- 
cus on WIMPs that dominantly undergo spin-dependent 
scattering, the limits on which are less stringent. 

We consider any model that can give rise to the fol- 
lowing effective axial interaction 
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between our DM particle, x, and standard model quarks, 
q, where Aq is a dimensionful parameter encoding the 
effective energy scale of the dynamics that generate this 
interaction. 

The X ^ P spin-dependent elastic scattering cross sec- 
tion at zero momentum transfer, resulting from the 
purely axial vector interaction, is given in [21j as 
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where the Aq factors are the spin fractions of the proton 
carried by a given quark [21H23J . It is this spin-dependent 
cross section that is central to our analysis. 

In addition to spin-dependent scattering, the DM 
model must have a suppressed annihilation cross section, 
compared to that required by successful DM genesis via 
thermal freeze-out, in order for DM particles to continu- 
ously accumulate in the Sun. The minimum suppression 
required for an effect to be seen is of the order of a p-wave 
suppression compared to the s-wave annihilation rate re- 
quired by standard freeze-out. Therefore, a suppression 
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quired in the standard small velocity expansion of the 
WIMP annihilation cross section, where v is the relative 



velocity between two colliding WIMP DM particles and 
the angled brackets represent a thermal average. 

A suppressed annihilation rate is not a generic feature 
of models of thermal freeze-out. However, models where 
the DM species possesses a particle-antiparticle asymme- 
try can lead to DM with zero annihilations today. The 
relic abundance in this case is assumed to be fixed by the 
value of the asymmetry, with the ratio of the baryon to 
DM relic abundance, flb/^x' determined by the dynam- 
ics that generate the asymmetry and is of order mp/m^. 
A number of attempts to link f2b to O^ have been made 
in the context of asymmetric DM (see e.g., [H]). Typ- 
ically, models of asymmetric DM involve DM particles 
with small masses and could therefore have implications 
for helioseismology. Previously, models of DM possessing 
large spin-dependent elastic scattering cross-section and 
an asymmetry have been investigated in the context of 
their effects on neutrino fluxes produced in the Sun |25j . 
Further related ideas can be found in [26" . 



WIMP CAPTURE IN THE SUN 

The accretion rate of the number of WIMPs captured 
by the Sun in the large spin-dependent scattering cross 
section limit, where all WIMPs intercepting the Sun are 
captured, is given by [27] 
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where fcsc. ~ 617kms~ is the escape velocity at the 
solar surface, v ~ 270 km s^^, ^(oo) ~ 0.75, and C = 1.77. 
This reduces to 
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assuming that the WIMP interaction cross section with 
matter is a^ > 10~^^ cm^ in order for all incoming 
WIMPs to undergo one or more scatterings while inside 
the Sun. We normalise the local DM density in the solar 
neighbourhood to pdm = 0.3GeVcm~^. Here, m-^ is the 
WIMP mass. Because of competing effects of annihila- 
tion and evaporation, the number of accreted WIMPs at 
time t is obtained by solving the differential equation 
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the self-annihilation 
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(fann.w) is tlic product of thcrmally-averaged WIMP self- 
annihilation cross section and velocity, and n^ = ^sm^ 
the number density of WIMPs inside the Sun, is as- 
sumed to be constant. The evaporation rate, Fovap., de- 
cays exponentially with temperature as ~ ^-GMm^/RT 
and is negligible with respect to the annihilation rate for 



m;^ > lOGeV. With this simpHfication, the population 
of WIMPs at time t is given by 



N{t) = (J"r) tanh(i/r)/2, 



(6) 



where the time-scale r = 1/\AFCa- For t ^ t, i.e. when 
the equilibrium between accretion and annihilation has 
been reached, the number of particles accreted is equal 
to the time-independent product Tt. 

Assuming that we are in the regime when the veloci- 
ties and positions of the WIMPs within the Sun follow a 
Maxwell-Boltzmann distribution with respect to its cen- 
tre, the amount of energy released by annihilation in a 
thermalisation volume centred in the compact star will 
have a radius EHl 
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Taking typical solar core conditions as having 
central internal temperature Tc = 1.57 x 10^ K 
and density /9c=154gcm~^, we find that Rth. — 

1/9 

0.03 (m^/lOGeV) ' Rq. We assume constant tem- 
perature and density within the WIMP thermalisation 
volume. Hence, WIMPs fill the solar core for masses 
in the range '-^S-lOGeV, and it is for this reason that 
there may be a significant imprint on helioseismology 
for both p- and g-modes that are sensitive to the sound 
crossing time across the solar interior. At lower masses, 
evaporation predominates, while at higher masses, the 
WIMPs within the Sun occupy a smaller volume. 



APPLICATION TO SOLAR MODELS 

The solar models shown here are evolved from the pre- 
main sequence using an updated version of the one- 
dimensional evolution codes described in Iben [29, 30, 31j. 
The evolution code uses the SIREFF EOS [sec^, Burg- 
ers (1969) [33] diffusion treatment as implemented by 
[34|, [35] . the nuclear reaction rates from Angulo et al. 
[36] with a correction to the ^*N rate from Formicola et 
al. [37] J and the OPAL opacities [3^ supplemented by 
the Ferguson et al. [39] low-temperature opacities. 

The models are calibrated to the present solar ra- 
dius Rq = 6.9599 X lO^^cm 00], luminosity Lq = 
3.846 X 10^3 ergs"^ [H], mass Mq = 1.989 x 10^3 g [42], 
and age 4.54 ± 0.04 Gyr [43]. Defining X and Z as the 
mass fraction of hydrogen, and the mass fraction of ele- 
ments heavier than helium, respectively, the models are 
calibrated to the photospheric ZjX ratio appropriate 



for either the Asplund, Grevesse and Sauval (2005) so- 
lar mixture [18] (hereafter AGS05), or the Grevesse and 
Noels (1993) solar mixture [44] (hereafter GN93). 

For the evolution models, the initial helium abundance, 
Yoj initial heavy element mass fraction, Zq, and mixing 
length to pressure-scale-height ratio, a are adjusted so 
that the final luminosity, radius, and surface ZjX match 
the above constraints to within uncertainties. 

From the final evolution model, a more finely zoned 
model is created for calculating the oscillation frequen- 
cies. The radial and non-radial non-adiabatic p-mode 
and g-mode frequencies are calculated using the La- 
grangian pulsation code developed by Pesnell [45] . (See 
[46] for additional references and description of the 
physics used in the evolution and pulsation codes and 
models.) 

The WIMP energy transport description is considered 
in two regimes, depending critically on the mean free 
path of the WIMPs and the scale radius of the system. 
This ratio is known as the Knudsen parameter, 
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where l^^i(r), the mean free path of the WIMPs relative 
to the rth element, and r^, the WIMP scale radius are 
defined, respectively, as 
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where <Ti is the elastic scattering cross section, Xi{r) is 
the mass fraction of isotopic species i at radius r (e.g. 
Xi is the proton), p{r) is the matter density in units of 
gcm~3 at radius r, and rrii is the mass of species i (e.g. 
the proton mass for hydrogen) and 
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where pc is the central solar density, T^ is the WIMP 
temperature, /cb is Boltzmann's constant, G is the grav- 
itational constant, and m^ is the WIMP mass. 

For the models considered here, only spin-dependent 
interactions are considered, with the contribution of hy- 
drogen overwhelming the spin interactions of the core. 
Thus, all formulas that are summed over the elements 
can be reduced to their hydrogen nuclei (i.e. proton) 
contribution. 

In [IT] HH] , a Monte Carlo method was used to solve 
several "generic toy star" models incorporating WIMPs 
in the conductive regime. The authors found that by 
using the conduction formula along with a suppression 
factor related to the Knudsen parameter, the entire non- 
local transport regime could be related to the conductive 
regime through the approximation 



Lx{r) = f{Kn)Lcond.{r) 



(11) 



where L^ is the total energy transferred from WIMPs 
to the nuclei in the intermediate regime between the 
conductive and non-local regimes, Lcond. is the energy 
transported by WIMPs in the conductive regime and the 
suppression factor /{K^) is given by 
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The Knudsen number Kq is the mean free path, in scale 
height units, that gives the most efhcient energy trans- 
port from the WIMPs within the Sun to the surrounding 
nuclei. Gould & Raffelt found this to be equal to ^ 0.4 
P71 H5] . For the solar models considered here, WIMPs 
are introduced into the energy transport by modifying 
the opacity 
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Here Krad.-i-e- is the combination of Rosseland mean ra- 
diative opacity and an effective opacity to take into ac- 
count electron thermal conduction that are added in re- 
ciprocal; Kcond. is the effective opacity derived by treating 
WIMP energy transport as a conductive process. 

In the introduction, we noted that for DM masses heav- 
ier than 5 — 10 GeV, direct detection places constraints on 
the spin-dependent elastic scattering cross section. How- 
ever, here we want to demonstrate and compare the ef- 
fects of our described DM particles on the Sun as we 
decrease their mass to values where their corresponding 
spin-dependent elastic scattering cross sections are per- 
mitted by direct detection limits. 

We have also explored models with WIMPs of differ- 
ent masses and interaction cross sections [[49 , and present 
this series as an illustration of the effects of WIMPs on 
a solar model when one parameter, namely, the WIMP 
mass, is varied. We chose an interaction cross section, 
o'x= 7x10"'^'^ cm^, and a very small annihilation cross 
section, (cTann.w) = 10^^" cm"^ s^^, to enhance the effect 
of WIMPS on the solar model. These results are in- 
tended to explore whether helioseismic signatures could 
have the potential to reveal or rule out the presence of 
a particular class of WIMP. At this stage we do not in- 
tend to provide rigorous helioseismic constraints on the 
properties of WIMPs. To this end, we compare the char- 
acteristic parameters of solar models based on the GN93 
and AGS05 solar abundances when including WIMPs of 
masses 50, 30, 20, 15, 10, and 5 GeV. 

The solar models including WIMPs use a tiny but non- 
zero annihilation rate. This is in contrast to the DM 
models outlined above where, due to the asymmetry in 
the DM species, the DM particles are unable to annihilate 
once captured by the Sun. This means that the annihila- 
tion rate used in our numerical solar models should also 
be zero. It turns out that the size of the effects mani- 
festing in solar properties plateaus such that the effect of 



decreasing the annihilation rate further does not signifi- 
cantly change the numerical results [1|. 



Effect on model structure 

In Table|l) we display the properties of standard solar 
models using either the GN93 or AGS05 abundances, 
as well as the properties of solar models when includ- 
ing WIMPs. Under the heading of "Model Calibration" 
we list the values of the following parameters: Xq, Yq 
and Zq are the initial mass fractions of H, He and metals 
(i.e. elements heavier than He); a is the mixing length- 
to-pressure scale height ratio; ZAMS is the zero-age main 
sequence; log (L/Lq) is the log luminosity in solar units; 
log (i?/i?0) is the log radius in solar units; Z/X (sur- 
face) is the surface ratio of metals to H mass fraction 
at the present solar age. Under "Solar Center Proper- 
ties" we list the values of the following parameters: T^. 
is the central temperature; pc is the central density; Kc 
is the central opacity; *B z/ flux is the predicted ®B neu- 
trino flux at Earth's surface, while the subsequent rows 
are the predicted total and ^B i/ fluxes for ■^^Cl detec- 
tors, in Solar Neutrino Units (SNUs), defined as 10~^^ 
absorptions per '^'^Cl atom per second. Under the head- 
ing "Helioseismology" , i?czB is the predicted ratio of the 
convection zone base to the solar radius, Rq, and Ycz is 
the predicted helium abundance in the convection zone. 
The constraints on these quantities from helioseismology 
are given in the table end-notes. 

Note that the model structure and calibration is con- 
siderably different for standard models without WIMPs 
calibrated to either the GN93 or AGS05 abundances. Us- 
ing the AGS05 abundances, which possess a smaller Z, a 
smaller helium mass fraction Y is required to compensate 
to increase the pressure in the core. Since more hydrogen 
fuel is available, both T^ and Pc are slightly reduced to 
produce the same luminosity. The location of the enve- 
lope convection base is determined by the radius where 
the temperature gradient exceeds the adiabatic gradient. 
For the AGS05 model, this point is reached at a lower 
temperature and larger radius because of the smaller frac- 
tion of heavier elements, particularly oxygen and neon, 
that are ionizing near the convection zone base and con- 
tributing to the opacity. 

We have added WIMPs to models calibrated to the 
AGS05 abundances. As we discuss below, WIMPs mainly 
would affect the innermost 10% of the Sun's radius that 
is sampled least well by the observed solar p-modes [5T] . 
For the models presented here, the effects become no- 
ticeable only for to^ < 20 GeV when the WIMPs are less 
tightly bound in their orbits around the solar centre and 
can transfer energy to larger distances from the core. 
With decreasing WIMP mass, WIMP energy transport 
cools and thermalizes the core to lower temperatures out 
to larger radii. Because the calibrated models need to 



TABLE I: Properties of Standard Solar Models and Solar Models including WIMPs" 



Model/WIMP Mass: 



GN93 AGS05 50 GeV 30 GeV 20 GeV 15 GeV 10 GeV 5 GeV 



Model Calibration:' 

Xo 

Yo 

Zo 

a 

Age-ZAMS (10® yrs) 

log (L/Lq) 

log (R/Re) 

Z/X (surface) 



Solar Center Properties at Solar Age: 

Tc (lO'^K) 

Pc (gcm-^) 

Y, 

Kc (cm^g-^) 

Sound Speedc (lO'^cms"^) 

**B i/flux (lO'^cm-^s-i) 

*B u flux ^^Cl detector (SNUs) 

Total u flux ^^Cl detector (SNUs) 



Helioseismology: 

RCZB {Rq) 

Ycz 



0.71000 


0.72950 


0.72950 


0.72950 


0.72950 


0.72950 


0.72950 


0.73570 


0.27027 


0.25698 


0.25698 


0.25698 


0.25693 


0.25690 


0.25677 


0.25057 


0.01973 


0.01352 


0.01352 


0.01352 


0.01357 


0.01360 


0.01373 


0.01373 


2.0423 


1.9916 


1.9913 


1.9910 


1.9963 


1.9990 


2.0121 


2.0734 


4.52 


4.52 


4.51 


4.50 


4.51 


4.51 


4.51 


4.51 


-5.74E-06 2.52E-06 


-9.31E-06 


7.86E-06 


-4.291E-06 


-7.16E-06 


-5.44E-06 


2.5E-05 


4.34E-07 


2.17E-06 


3.04E-06 


-1.74E-06 


3.04E-06 


8.69E-07 


2.61E-06 


3.04E-06 


0.0246 


0.01628 


0.01628 


0.01629 


0.01635 


0.01639 


0.01657 


0.01654 


15.64 


15.42 


15.172 


15.005 


14.814 


14.646 


14.330 


13.482 


152.40 


148.96 


149.84 


150.70 


152.45 


154.11 


158.31 


175.57 


0.6329 


0.6183 


0.6076 


0.6009 


0.5930 


0.5851 


0.5780 


0.5488 


1.231 


1.261 


0.02408 


0.02254 


0.02296 


0.02290 


0.02183 


0.01756 


5.083 


5.058 


5.071 


5.061 


5.050 


5.040 


5.007 


4.94 


5.26 


4.30 


4.19 


3.98 


3.68 


3.30 


2.60 


1.04 


5.99 


4.91 


4.77 


4.53 


4.19 


3.76 


2.96 


1.18 


7.60 


6.32 


6.18 


5.92 


5.56 


5.10 


4.23 


2.21 


0.7133 


0.7294 


0.7294 


0.7293 


0.7280 


0.7280 


0.7275 


0.7220 


0.2419 


0.2273 


0.2273 


0.2273 


0.2273 


0.2273 


0.2274 


0.2227 



" Characteristic parameters of solar models based on the GN93 and AGS05 abundances, and AGS05 models including WIMPs of masses 
50, 30, 20, 15, 10, and 5 GeV. All models were run with (cann'') = 1 x 10"'*'' cm^ s~^ and a spin-dependent interaction cross section of 
(Ta; = 7 X 10~^^ cm^ in order to enhanced the effects WIMPs have on the solar model. 

''Xo, loi aud Zq are the initial mass fractions of hydrogen, helium and elements heavier than H and He, respectively; a is the mixing 
length to pressure-scale-height ratio; ZAMS is the zero-age main sequence. 

'^Tc, Pc Yc, Kc, and Sound Speedc are the central temperature, density, helium mass fraction, opacity, and sound speed, respectively; 
the neutrino fluxes are given for the *B and total fluxes at Earth's surface, in cin~^s~^ and in Solar Neutrino Units (SNUs) of 10~^® 
absorptions per '^''Cl atom per second. 

''The seismically-inferred CZ helium mass fraction and CZ base radius are 0.248 it 0.003 and 0.713 it 0.001 Rq, respectively |50j . 



remain in hydrostatic equilibrium and to generate the 
same luminosity, pc and the central hydrogen abundance 
increase to compensate for the cooler temperatures. The 
lower temperatures reduce the predicted neutrino flux, 
particularly the *B v flux, which has a steep tempera- 
ture dependence of T^^ near the solar center [55^ . While 
the current solar neutrino experiments can accommodate 
the ^B neutrino flux predicted by standard solar models 
using either the GN93 or AGS05 abundances, more work 
is required to determine whether they can accommodate 
a flux as low as predicted for the discussed solar models 
including WIMPs. 

With increasing density concentration in the solar core, 
the envelope becomes less condensed, and would nor- 
mally have a larger radius. Therefore, the mixing length 
a is slightly increased to calibrate the model to the ob- 



served solar radius, resulting in increased convective effi- 
ciency and the onset of convection occurring at a slightly 
smaller radius. This change in a is nevertheless much too 
small to deepen the convection zone base radius -RczBj 
while retaining the AGS05 abundances, to the value of 
0.713 ± 0.001 Rq determined from helioseismic data. 

In Fig.jT] we plot the solar temperature as a function of 
fractional solar radius, for solar models with WIMPs with 
masses vn^ = 5 — 50 GeV, compared with the tempera- 
ture profiles corresponding to the AGS05 and GN93 mod- 
els. As discussed above, the average orbital radius and in- 
teraction region of the WIMPs increases with decreasing 
mass. As seen in Table|lj the WIMPs reduce the effective 
opacity in the core to only approximately 0.02cm^g~^, 
compared to 1.2cm^g~^ for the standard solar models. 
The transport of energy by the WIMPs is so efficient that 
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FIG. 1: Central temperature versus fractional solar radius for 
solar models with WIMPS with masses 5-50 GeV compared to 
the AGS05 and GN93 model temperature profile. All WIMP 
models were run with ((Tann.w) = 10"*" cm'^ (annihilations 

-- 7 X 



suppressed) and a spin-dependent cross section a 
10"''^ cm^, using the AGS05 element abundances. 



the core essentially becomes isothermal out to the edge 
of the interaction region, where the temperature gradi- 
ent approaches that of the standard model when WIMPs 
are omitted. For the most extreme 5 GeV WIMP mass 
model, the solar temperature is significantly reduced out 
to a radius of approximately 0.1 Rq. 



Effect on sound speed 

We now address the question as to whether the structural 
changes in the solar core discussed in the previous section 
lead to a detectable helioseismic signature. Fig. [2] shows 
the difference between the sound speed inferred from so- 
lar p-modes, compared to that generated by standard 
models as well as those models including WIMPs with 
5 < TO^ < 20 GeV. 

Because only a few solar p-mode eigenfunctions of the 
lowest degree < £ < 2 have significant amplitude near 
the solar core |51) . the sound-speed inversions using p- 
modes are not sensitive to solar central conditions with 
high accuracy for radii within 0.06 i?0 [S3]. In Fig.[2J we 
omit plotting the results corresponding to m-^ = 30 and 
50 GeV, as they nearly coincide with those from the stan- 
dard model with AGS05 abundances. These sound-speed 
difference curves nearly coincide because such WIMPs 
only affect the model structure and sound speed profile 
for radii within 0.04 i?0. 

From Fig.[2J we observed that both 10 GeVand 
20GeVWIMPs have only a small effect on the sound- 
speed profile outside the central core at radii 0.06 < r < 
0.2 _R0. The models with WIMPs predict a lower cen- 



FIG. 2: Relative difference between the inferred and cal- 
culated sound speeds for solar models constructed with the 
GN93 and AGS05 abundances. Also compared are models 
with energy transport by WIMPs, based on the AGS05 abun- 
dances, for masses rUx = 5,10,15 and 20 GeV, (tJann.u) = 
10"^" cm'' s~^ and a spin-dependent interaction cross section 
oi ax = 7 X 10"'^^ cm^. The inferred sound speed is from |53) . 
This combination of parameters enhances the effect WIMPs 
have on the sound-speed profile. 



tral helium abundance (see Table^ and a correspond- 
ing lower central mean molecular weight, /i, while at the 
same time the temperature profile becomes identical to 
that of the standard model at radius r > 0.08 i?0 for 
m^ = 10 GeV, and > 0.06 Rq for m^ = 20 GeV. There- 
fore, the sound speed, which is proportional to y/{T/ fi), 
is increased in this region for the WIMP models. Un- 
fortunately, the small increase in sound speed in this re- 
gion, that could be diagnosed by sound speed inversions, 
is in the wrong direction to reduce the discrepancy with 
the seismic inversions observed for the standard model 
AGS05 abundances. 

We also observe from Fig.|2J that we can use the 
discrepancy in predictions for core sound-speed and p- 
mode frequencies of the solar model including 5 GeVmass 
WIMPs to rule out this model. The large core tem- 
perature decrease produced by the WIMPs extends to 
10% of the solar radius, far enough out for p-mode 
sensitivity. The lower temperature results in a slower 
sound speed compared to that predicted by the standard 
AGS05 model for r < 0.1 i?0. When the temperature 
profile joins that of the AGS05 model at 0.1 i?0, the 
much-reduced mean molecular weight results in an in- 
crease in sound speed compared to the standard AGS05 
model at larger radii out to 0.4 _R0. Between 0.4 i?0 and 
RczB = 0.722 i?0, the 5 GeV WIMP mass model actu- 
ally mitigates the sound speed discrepancy for the AGS05 
abundances, because the lower Y abundance of the cali- 
brated model increases the sound speed, and the slightly 
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FIG. 3: Difference between the calculated and observed small 
frequency separations for p-modes of degree £ — and £ = 
2 against radial order for 1 = modes. The observations 
are from [51], corrected for solar-cycle variations using the 
method described in 1551. 



larger mixing length/pressure-scale-height ratio improves 
the efficiency of convection and slightly deepens the con- 
vection zone. However, for this model and all of the 
AGS05 abundance models presented here with/without 
WIMPs, the Y abundance needed to calibrate the model 
to the solar luminosity is too low compared to that de- 
rived from the signature of helium ionization in the con- 
vection zone [50] , and the convection zones of these mod- 
els are also still too shallow compared to the seismically- 
derived value. 



Effect on small frequency separations 

The differences between the very precisely measured £ = 
and i = 2 p-mode frequencies offset by one radial or- 
der, known as the small separations, are somewhat more 
sensitive to the solar central conditions. 

Fig. [3] shows the difference between the calculated and 
observed small separations. Notice that the calculated 
small separations for the standard model with GN93 
abundances agree well with observations, whereas the 
separations for the model with the AGS05 abundances 
are about 0.2 fjHz too high. There is very little effect 
on these separations for the AGSG5 abundance models 
including WIMPs with mass m^ > 10 GeV. The effect 
becomes more pronounced for the 10 GeV WIMP mass 
model, where the WIMPs introduce a decreasing slope 
with increasing radial order. This slope becomes much 
more extreme for the 5 GeV model, and consequently we 
find that helioseismic observations would rule out such a 



model. 



Effect on gravity modes 

While the sound speed inferred from solar pressure modes 
and other solar structure properties outside the solar core 
are not very sensitive to the effects of WIMPs, the longer- 
period solar gravity modes, which have largest ampli- 
tudes in the solar core but very small amplitude near the 
solar surface [5T], show more sensitivity. To date, only 
one or at most a few g-modes have been possibly de- 
tected in long-term observations by the SOHO spacecraft 
[551 [57] . Tableln] lists the calculated g-mode frequencies 
for the standard (i.e., without WIMPs) solar models with 
GN93 and AGS05 abundances, and the calculated differ- 
ences between the frequencies of the AGS05 model and 
the frequencies of the models with WIMP masses of 20, 
10, and 5 GeV. The g-modes with the lowest-degree i 
have the largest predicted amplitudes at the solar surface, 
and therefore are the most likely to be detected. As can 
be observed from Fig.lSJ the changes in frequency when 
WIMPs of these masses are included are several to tens 
of fjUz. If g-mode frequencies are isolated and the modes 
are identified unambiguously, the frequencies should be 
measurable to less than a /xHz, and so will easily allow 
one to distinguish between solar models with/without 
WIMPs, or between WIMP models with different values 

of TOv. 



DISCUSSION AND CONCLUSIONS 

Our main conclusion is that both direct detection and 
accelerator probes may be complemented by using the 
Sun as a probe of DM. Models of DM that have large 
spin-dependent interactions and an intrinsic asymmetry 
that prevents post freeze-out annihilations can signifi- 
cantly lower the central temperature of the Sun as well 
as the resulting ^B neutrino flux. For WIMP masses 
m-^ > 10 GeV, the presence of WIMPs does not signifi- 
cantly affect currently available helioseismic constraints. 
However, for WIMP masses of 10 GeV or lighter, con- 
straints on sound speed and small frequency separations 
between i = and £ = 2 p-modes can begin to constrain 
and rule out the presence of WIMPs with the cross sec- 
tions utilised here. 

Our study is motivated in a large part by the re- 
cently revised solar abundances [THl [IS] which result in 
solar models that cannot reproduce the currently ob- 
served helioseismic data, with numerous attempts to re- 
store agreement being met with only partial success (see 
e.g., [SH1[53). This means that additional physics must 
be incorporated into solar modelling, and dark matter is 
among the options that merit detailed consideration. 



TABLE II: Comparison of GN93, AGS05 and WIMP Model g-Mode Predictions' 
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Typical low-order g-mode predictions for the GN93 and AGS05 models, and predicted difference from the AGS05 model prediction for 



solar models with WIMPs of va^ = 20, 10, and 5 GeV. All WIMP models were run with ((Tann'u) 
cm^ in order to enhance the effects on the solar model. The units of u and the quantity Au"^^ 
are fjUz. 
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Since the original submission of our paper in May 2010, 
an additional paper appeared on solar models including 
WIMPs and the implications for reconciling the new so- 
lar abundances with helioscismology [60j . In agreement 
with [60] , our explorations to date do not show any real- 
istic path in which the inclusion of WIMPs will mitigate 
this problem. Even for the large interaction cross sec- 
tion and small annihilation cross section considered here, 
the inclusion of WIMPs of mass m^ > 10 GeV has lit- 
tle effect on presently observable helioseismic signatures. 
The inclusion of WIMPs with masses of 10 GeV or lighter 
worsens the agreement with the helioseismically-inferred 
sound speed at radii 0.1 < r < 0.2 Rq, only shghtly deep- 
ens the predicted convection-zone depth, and introduces 
a trend with radial order in the low-degree p-mode small 
separations that is not observed in the data. Our pri- 
mary new result is that WIMP masses of ~ 5 GeV may 
be excluded for spin-dependent interactions in a specified 
cross-section range, thereby complementing direct detec- 
tion experiments in a region that they access only with 
great difhculty provided the WIMPs annihilation cross 
section is suppressed. 

While here we do not discuss whether these WIMP 
models could accommodate measurements of the *B neu- 
trino flux from current solar neutrino experiments such 
as Super-Kamiokande III [61] , SNO [62] or Borexino [63] , 
with precisions of ^10% and theoretical expectations of 
up to ^--^20% depending on the solar composition [I] |5D] , 
a more detailed study of the low-mass region of WIMP 



parameter space and its consistency with current exper- 
imental data is deferred to a later paper. There we will 
address the question of whether future helioseismic ob- 
servations, most notably using g-modes, and solar neu- 
trinos, may be able to constrain the allowable DM pa- 
rameter space in a mass range that is of current interest 
for direct detection. 

Finally, we note that for solar mass stars near the cen- 
tre of the Galaxy, where the WIMP density is enhanced 
by up to some 6 orders of magnitude relative to that in 
the solar neighbourhood, the effect of the redistribution 
of energy in the stellar core may generate a significant 
reduction of the main-sequence lifetimes. We leave an 
investigation of this scenario to our future work. 
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